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We find that diverse nonlinear waves such as soliton, Akhmediev breather, and rogue waves, can 
emerge and interplay with each other in a two-mode coupled system. It provides a good platform 
to study interaction between different kinds of nonlinear waves. In particular, we obtain dark rogue 
waves analytically for the first time in the coupled system, and find that two rogue waves can appear 
in the temporal-spatial distribution. Possible ways to observe these nonlinear waves are discussed. 
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I. INTRODUCTION 

Solitons are localized waves arising from the interplay 
between self-focusing (self-defocusing) and dispersion ef¬ 
fect, and could be one of the intense studies in nonlinear 
science. The well-known solitons of the scalar nonlin¬ 
ear Schrddinger equation (NLSE) have been studied in 
many different systems including plasmas, optical fibers 
and cold atoms, mainly including bright soltons 
and dark solitons |5l-[lQj. In addition to solitons, the 
NLSE admits other classes of localized structure, such as 
Akhmediev breather (AB) 0, and rogue waves (RW) 
0 . Moreover, the AB and RW have been observed in 
one-mode nonlinear fibers experimentally 0. 

Recently, the studies for solitons have been extended 
to multi-component coupled systems. It has been 
reported that the bright-bright(B-B), dark-dark(D-D), 
bright-dark (B-^ solitons can exist in different param¬ 
eters regimes |l3-[lq. The bright-dark soliton dynamics 
have been observed experimentally in [T^ . Besides these 
solitons, studies have indicated that there exist other 
kinds of nonlinear localized waves in the coupled system, 
such as AB 00 , and vector RWs 0 . These stud¬ 
ies indicate that nonlinear waves in coupled system are 
much more diverse than the ones in uncoupled systems. 

In this paper, we study diverse families of nonlinear 
localized waves in a two-mode optical fiber. The proper 
conditions for these nonlinear waves are given explicitly. 
Notably, interactions between different kinds of nonlin¬ 
ear waves can be observed analytically. Eor example, 
bright or dark soliton interplay with RW could be ob¬ 
served, which provides a good platform to study inter¬ 
action between RW and other nonlinear waves. More¬ 
over, we prove that dark RW indeed can exist in one 
component of the coupled system, and report that two 
RWs could appear in temporal-spatial distribution. Es¬ 
pecially, there are some additional requirements on the 
nonzero backgrounds for vector RWs, which is quite dif¬ 
ferent from the RW in single-component nonlinear sys- 
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tern. As an example to demonstrate possibilities to ob¬ 
serve these nonlinear waves, we present one possible way 
to create RW in the two-mode nonlinear fiber system. 

The paper is organized as follows. In Section II, we 
present three families of nonlinear localized wave solu¬ 
tions and discuss the explicit conditions under which they 
could exist in detail. In Section III, the possibilities to 
observe them are discussed. The conclusion and discus¬ 
sion are made in Section IV. 


II. ANALYTICAL VECTOR NONLINEAR 
WAVE SOLUTIONS 

It is well known that coupled NLS equations are often 
used to describe the interaction among the modes in non¬ 
linear optics. We begin with the well known two-coupled 
NLSE in dimensionless form 

^2 + 2 ^ 2 |' 02|^]'01 = 0 , ( 1 ) 

+ ^2 ^2 + + 2 ^ 2 |' 02|^]'02 = 0 , ( 2 ) 

where and -02 are complex envelopes of the electric 
fields of the two modes, z denotes propagation distance, 
and t represents the retarded time. The ai and (72 depend 
on the signs of the group velocity dispersion(GVD) in 
each mode, i.e., a = +1 for anomalous GVD and a = —1 
for normal GVD. Here, we just consider the anomalous 
GVD case, ai = (72 = 1. ^ 1,^2 are nonlinear parame¬ 
ters determined by properties of the Kerr medium with 
electrostrictive mechanism of nonlinearity [ 2 ^. When 
9i = 92 = Q, it will become the well-known Manakov 
model [3: • 

The Eq.(l) has been solved to get soliton solution on 
trivial background through Hirota bilinear method in 
0 . Performing Darboux-transformation from a triv¬ 
ial seed solution, one could get the bright-bright soli¬ 
tons [TsI . It has been reported that solitons could collide 
inelastically and there are shape-changing collisions for 
coupled system, which are different from uncoupled sys¬ 
tem 0 . However, it is not possible to study B-D soliton, 
ABs or RWs on trivial background. Next, we will solve 
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it from nontrivial seed solutions. The nontrivial seed so¬ 
lutions are derived as follows 


V^io 

^^20 


Si exp [iOi[t, z]] 

51 exp [ikit + i{2gis‘l + 2g2sl — ikl)z], (3) 

5 2 exp [i 02 [t,z]] 

S 2 exp [ik 2 t + i{2gisl + 2 ^ 2 «S 2 - ikl)z], (4) 


where Si and S 2 are two arbitrary real constants, and de¬ 
note the backgrounds in which localized nonlinear waves 
emerge. ki and k 2 denote the frequencies of the plane 
wave background in the two modes respectively. We 
will solve the coupled system analytically through Dar- 
boux transformation method, which can transfer nonlin¬ 
ear problem to linear one. The corresponding Lax-pair 
of Eq.(l) with ^ 1,^2 > 0 can be derived as 



where 

/ V^V’2 \ 

u = |A 0 , 

V-V^^2 0 iA J 


v = ux + 


-iy'gig2'lp2^1 


iy/^'<p2t _ 

-^^/gmi’2^l 

J3 


and 


Jl = +*5'2|^/’2p 

J 2 = 

J 3 = -i92\'>p2\'^ ■ 


Hereafter, the overbar denotes the complex conjugate. 
With spectral parameter A and the nontrivial seed solu¬ 
tions, one can solve the following equation of r to get the 
diagonalized or Jordan forms of U and V through the 
method presented in H- Then, one could get nonlinear 
wave solutions through performing Darboux transforma¬ 
tion method. The equation of r, which is the eigenvalue 
equation of the transformed [/, is calculated as 

+ aiT + 61 = 0, (7) 

where 

ai = be — ab — ac ^ giSi ^ g2S2^ 
bi = abc- gislc- g2slb, 

= '<^2—+-(/ci +/C 2 ), 

6 = ~ ^ 2 ), 

c = i^ + |(2fc2 - ki). 


There will be three cases for the roots rj. Each one 
corresponds to a different family of nonlinear waves in 
the coupled systems. 

Casel: The three roots are all different, we find that 
there are bright-dark solitons pair interplay with Akhme- 
diev breather (B-DAB), and Akhmediev breathers- 
Akhmediev breathers (ABs-ABs) solutions in the coupled 
system. 

Case 2 : There are one single root and double roots, 
the bright-dark interplay with rogue waves (B-DRW), 
Akhmediev breather-Akhmediev breather interplay with 
RW (AB-ABRW) solutions could exist in the coupled sys¬ 
tem. 

Case3: There are triple roots, the RWs with no other 
type nonlinear waves solutions could exist with certain 
conditions. 

We find that there are many explicit requirements on 
signals or backgrounds for different kinds of vector non¬ 
linear waves. It should be pointed out that similar stud¬ 
ies have been done on nontrivial backgrounds in [ 13 , I 2 Q| - 
[ 22 [ . Distinct from them, we investigate all possible cases 
and report explicit conditions for these nonlinear wave 
solutions in the nonlinear coupled systems. The whole 
analysis would help us to understand the relations be¬ 
tween vector solitons, breathers, and RWs. Eurt her more, 
we report that there are some new types of nonlinear 
waves, such as dark rouge waves, just two rouge waves in 
temporal-spatial distribution. Next, we will discuss these 
three circumstances. 


A. Bright-dark solitons interplay with Akhmediev 
breather, and Akhmediev breathers-Akhmediev 
breathers 


When Tj (j=l,2,3) are three different roots for Eq.(7), 
we can solve the Lax-pair to get 4>i, 4 > 2 , 4>3 as 


^l[t,z] = 


(^>01 + ^02 + ^ 03 ) 
I 


xExp 


^2[t,z] = - 


3 


{0i[t,z] + 92[t,z]) 


(8) 


xExp 


Ti — b 

I 


$01-r$02-r^03 


^3[t,z] = 


3 

V^S2 


T2 -b Ts-b 

{92[t,z] - 29i[t,z]) 

V 92 S 2 


(9) 




Ti - C 


01 


T2 - C 


^ y^S2 
^02 -^03 


T3 - C 


xExp 


-i9i[t,z] - 292[t,z]) 


( 10 ) 
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where 

^oi[t,z] = AiEx.p[rit ^ Ir^z2{X - ki - k2)riz/3] 
xExp[F[z]], 

^02[t,z] = A2Ex.p[r2tFlr2Z-\-2{X-ki-k2)r2z/3] 
xExp[F[z]], 

^03[t,z] = A3Exp[r3tFlrizF2{X-ki-k2)T3z/3] 
xExp[F[z]], 

where A = aoFibo and F[z] = 21 /9{X‘^F {ki Fk 2 )XFkf F 
k^ — kik 2 F 3gis‘l F 3g2S2)z. Between these expressions, 
the parameters ao, bo, Aj{j = 1,2,3), and si, S 2 are all 
real numbers which relate with the initial condition of 
soliton or AB, such as initial coordinate , initial velocity, 
and initial shape. Then, one can perform the following 
Darboux transformation to get one family of solutions for 
Eq.(l) and (2) 


V’l 

— — 

1 i{X — X)u 

( 11 ) 

V^l+|w|2 + |t;|2’ 

V’2 

= i^20 — 

1 i{X — X)v 

( 12 ) 

V^l+|w|2 + |t;|2’ 


where u = ^ and 'T' = ^. This is a generic solution 
which can be used to study the properties of B-D solitons, 
AB-AB, B-AB, and D-AB for coupled nonlinear systems. 
To get B-B solitons, it is much easier to solve the Lax-pair 
form trivial seed solution, namely, ^lo = i ^20 = 0. Based 
on the generic solution, it is more convenient to study 
their dynamics on different nontrivial backgrounds for 
the coupled systems. We find that there are much more 
abundant nonlinear waves on the nontrivial background 
than the ones on the trivial background. 



FIG. 1: (color online) The evolution plot of B-DAB in the 
coupled system, (a) for bright soliton in ^|;l component and 
(b) for dark with an AB in 'b ’2 component. It is seen that 
the bright soliton in one field is reflected like hard wall by 
the AB in the other one. The dark soliton is reflected well by 
the AB too. The coefficients are ao = 0, 5o = 1-2, = ^2 = 

0.25, Si = 0.001, S2 = 2, Ai = 1,^2 = 4,^3 = 3,/ci = —0.5, 
and k 2 = 0.5. 

When ^ 0, 52 7 ^ 0, with the certain parameters, 
the bight-dark interplay with AB (B-DAB) could exist, 
such as Fig. 1. Particularly, Fig. 1(a) shows that one 
bright soliton in 'ipi component is reflected by the other 


one’s density distribution. Considering that the density 
distribution of one component could be seen as nonlin¬ 
ear potential distribution for solitons in the other field, 
it is not hard to understand the refection effect. Corre¬ 
spondingly, one can observe a dark soliton collide with a 
AB and its shape is changed in Fig. 1(b). The reflect¬ 
ing inelastic collisions happen between nonlinear waves 
with different structures. Considering that Becker et al. 
have observed dark soliton collide with the B-D soliton 
in BEC [a, it is believed there are some possibilities to 
observe the refection collision presented here. 



FIG. 2: (color online) The evolution plot of ABs-ABs in two 
components coupled system, (a) for ABs soliton in '0i com¬ 
ponent and (b) for ABs in '02 component. It is shown that 
one AB collides with a AB inelastically. The coefficients are 
ao = 0, 5o = 1.2, gi — g 2 — 0.25, Si = 1, S 2 = 1, Ai = 1, A 2 = 
4, A 3 = 3, /ci = —0.5, and k 2 = 0.5. 

When 5 i 7 ^ 0 and 52 7 ^ 0, there are many kinds of dif¬ 
ferent soliton and AB appearing in the coupled system. 
In the nontrivial background, ABs-ABs could exist. The 
reflection between AB and AB can be observed, just as 
shown in Fig. 2(a). Two ABs could collide inelastically 
and then merge into one AB, as shown in Fig. 2(b). In¬ 
elastic collision for dark, bright, AB with AB could exist 
under some certain conditions. These characters could 
be used to design optical switches in nonlinear optics. 


B. Bright-dark solitons or Akhmediev Breather 
interplay with Rogue Wave 

It is found that when the equation of r has double 
roots, namely, ri = — 2 t 2 and T 2 = T 3 , there is a certain 
requirement on the spectral parameter A. The parame¬ 
ter A which mainly determines the form of initial signals 
should satisfy the equation as follows 

[3{ki F ^ 2 )^ + 4((7 — A)]X'^ 

+ [4(A:i + k2)^ F {ki F k2){SC — 6 A) — 45]A^ 

+ [4(A^i + k2)‘^C F 4G^/3 — A^ — 6{ki F k2)B]X^ 

F[4:{ki F k2)C‘^l3 — 2AB]X H——- B^ = 0, (13) 
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where 

A = 2(2ki — k2)(2k2 — ki) + (ki + k2)^ + 

B = 9giSi{2k2 — ki) + {ki + k2){2ki — k2){2k2 — ki) 
+9g2sl{2ki - k2), 

C = {ki + ^ 2 )^ ~ i^ki — k2){2k2 — ki) + 9g\s\ + 9^2S2- 

Under this requirement condition, the roots of r will be 
given as 


Tl 


- 2 r 2 , 


T2 = ts = 


Hi{\) 

H2{\y 


(14) 

(15) 


where 


Hx{X) = I[2X^+2.{kx + k2)X^+ {2{2kx-k2){2k2-kx) 
+(fci + ^ 2 )^ + 9gis\ + 9g2s\)X 
+9g2sl{2ki - k2) + 9gis\{2k2 - ki) 

+{ki + k2)(2ki - k2)(2k2 - ki)], 

H2{X) = 6 A^ + 6 (fci + k2)X + 2 (fci + ^ 2 )^ 

-2{2ki - k2){2k2 - ki) + I 851 S 1 + 185252 - 


RW emerges in the second component, which can be well 
shown by density distribution plot near z = 9. Moreover, 
the relative location between bright or dark solitons and 
RW can be varied through changing the parameters Aj. 
Especially, when = 0, the AB-ABRW solution will 
correspond to pure one vector RW solution, shown in Fig. 
4(c) and (d). In Fig. 4(d), it is shown well that the dark 
RW indeed exists in the second component. The simi¬ 
lar dark RWs have been expected in mixed BEC system 
through numeric stimulation in [^. Namely, we could 
obtain analytical dark RW solution with the correspond¬ 
ing parameters in Fig. 4(d). Furthermore, when A 3 = 0, 
the solution will become AB-AB solution with no RW 
any more. The parameter A 2 just affects the place where 
RW appears. 





FIG. 3: (color online) The evolution of B-DRW in the cou¬ 
pled system, (a) for one bright soliton with RW in com¬ 
ponent and (b) for one dark with a RW in '02 component. It 
is shown that one dark soliton collides with RW elastically. 
Bright soliton are attracted by RW when being next to it. 
The coefficients are ao = —0.1,5o = —2,^i = ^2 = 0.25, si = 
0.01, S 2 = 2, Ai = 1, A 2 = 4, As = 3, ki — —1, and k 2 = 0.1. 



FIG. 4: (color online) The plot of AB-ABRW in the coupled 
system, (a) for AB with a RW in 0i component and (b) for 
AB with a dark RW in 02 component. For (a) and (b), Ai = 
1, A 2 =4. (c) The plot of one RW in 0i component and (d) 
The plot of one RW in 02 component. It is seen that the RW 
in 02 component is a dark one. (c) and (d) with Ai = A 2 = 0 . 
The other coefficients are ao = —0.53, bo = 2.26, gi — g 2 — 
0.5, Si = 1, S 2 = 2 , As = 3, /ci = —1, and k 2 = 0.1. 


Under these conditions, one RW appears, namely, the 
AB in previous part will evolve as a RW. And the RW 
with other type nonlinear waves solution is given in Ap¬ 
pendix A. One can see that the form of solution partly 
includes rational functions. When the amplitude of one 
component’s background Si ^ 0 , and S 2 > 0 , the bright- 
dark interplay with RW (B-DRW) could be observed, 
such as Fig. 3. It is seen from Fig. 3(a) that the bright 
soliton is attracted by the RW, and collide with it, then 
the shape of dark soliton is unchanged after RW disap¬ 
pearing. When si,S 2 > 0, the ABs emerge with RW 
(AB-ABRW) in each component, such as Fig. 4. One 
normal RW appears around t = 0 place near 2 : = 0 in 
the first component field, and correspondingly, one dark 


On arbitrary nontrivial background, when the initial 
condition can be given as nearly with the certain A be¬ 
fore, there are many possibilities to observe one RW in¬ 
teracted with bright, dark, or ABs. These interaction 
phenomena are impossible to be observed for scalar NFS 
system, for bright, dark solitons and RW can not exist si¬ 
multaneously with the same conditions. The generalized 
nonlinear waves solution provides us a good platform to 
study the interactions between RW and other type non¬ 
linear waves. On the other hand, the nonlinear wave 
solution partly including rational functions corresponds 
to a RW with other type nonlinear waves. When Ai = 0, 
the AB-ABRW solution will be all rational form and it is 
one vector RW solution. One could guess that the vector 
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RWs with no other type nonlinear waves solution could 
have the all rational form, just like RW solution in uncou¬ 
pled system. Then, we will try to derive the all rational 
form solutions. 


C. Rogue Waves with no other nonlinear waves 

When r has triple roots, we can derive the solutions 
with all rational functions, which usually correspond to 
RWs phenomena [HI Interestingly, we find that 

there are some certain requirements on the nontrivial 
backgrounds of the two components for these RWs so¬ 
lutions. The conditions are required as 


gisj = 92 sI, 


(16) 

(17) 


\ki -k2\ = x/^si, 
which means that there are requirement on the non- 



FIG. 5: (color online) The evolution plot of only one RW 
in coupled system, (a) for one RW in component, (b) for 
one RW in ^2 component. The coefficients are = ^2 = 
0.25, Si = S 2 = l,Ai = 1,^2 = 1,^3 = 0,ki = 0.55, and 
k2 = 0.05. 

linear parameter and the amplitude of each component, 
and the difference of their wave vectors should be related 
with nonlinear parameter and the amplitude in a certain 
way. Then, RWs without any other type waves, could 
be observed possibly on the certain backgrounds. Mean¬ 
while, the ideal initial conditions for the RW are pre¬ 
sented in the coupled system. The parameter A which 
mainly determines the initial shape of vector nonlinear 
waves, should be set with 

A--^-(18) 

With these certain conditions, the generic form of vec¬ 
tor RWs could be given as 

i + | y |2 + |^|2 > uy; 

_ n I ^V^siW2it,z) 

^2 - [ 1 + 1 + H 2 + |^|2 ^ ( 20 ) 

where u = ITi(t, 2:)y^Si, and v = W 2 {t, z)y^S 2 - 
The expressions of Wi^ 2 {t,z) are all rational functions. 



FIG. 6: (color online) The evolution plot of two RWs in cou¬ 
pled system, (a) for two similar RWs in component, (b) 
for two similar RWs in ^2 component. The coefficients are 
gi = g 2 = 0.25, Si = S 2 = 4, = 1,^2 = 1,^3 = 0.2,= 

2.05, and k 2 = 0.05. 


presented in Appendix B. Therefore, the vector waves 
solution could be vector RWs solution, which can be ver¬ 
ified by the following RWs plots. When gi = g 2 = 
the RWs solution could become the RWs in [^. When 
A3 = 0, only one RW can be observed in both compo¬ 
nents, shown in Fig. 5. The density of them are distinct 
in temporal-spatial distribution, and their structures are 
similar to the well-known RW in single-component sys¬ 
tem. When A3 7^ 0, there are two RWs appearing in the 
temporal-spatial distribution, shown in Fig. 6. There 
are two RWs appearing in temporal-spatial distribution, 
which is very distinct from the higher-order RW in un¬ 
coupled system. In the uncoupled systems, it is nor pos¬ 
sible to observe just two RWs appearing in the whole 
temporal-spatial distribution even for higher-order RWs 

[21,123. 

When parameters Aj vary, namely, the initial nonlin¬ 
ear waves are slightly different, the RWs properties will 
be changed. For example, when they are next to each 
other, the classical shape of RW [^,[^, called eye shape, 
could be changed, such as fig. 7. However, the RW solu¬ 
tion form is an ideal condition. It could be very hard to 
realize them precisely. From the generic solution Eq.(ll) 
and (12), when the initial condition is a bit different from 
the ideal initial condition, one can observe the character 
of their evolution, such as Fig. 8. Comparing Fig. 6 and 
Fig. 8, we can know that the vector RW could be seen 
as a limit case of the general solution under the certain 
requirement condition. When the initial conditions are 
similar to the ideal one, the evolution of them will be 
close to the vector RWs. This is similar to the case that 
RW is the limit of AB in the uncoupled system [T^ . 


III. POSSIBILITIES TO OBSERVE THESE 
VECTOR NONLINEAR WAVES 

Considering the experiments on RW in nonlinear fibers 
with anomalous GVD [T^, [^, , which have shown that 

the simple scalar NLS could describe nonlinear waves in 
nonlinear fibers well, we expect that these different vector 









6 



FIG. 7: (color online) The evolution plot of two vector RWs 
in the coupled system, (a) for two different RWs in '0i com¬ 
ponent, (b) for two different RWs in component. It is seen 
that dark RW could be observed in one component of the 
coupled system. The coefficients are gi — g 2 — 0.5, si = S 2 = 
2, Ai = 0,A2 = 0, As = 158, ki = 0.807, and k 2 = 0.05. 



FIG. 8: (Golor online) The evolution of whole ABs-ABs 
density in the coupled system. The initial condition ap¬ 
proaches the ideal condition for RWs, and the corresponding 
required backgrounds are set. It is seen that the evolution of 
them is close to RWs in Fig.6. The explicit coefficients are 
ao = —1.05, 5o = 5.196,5'! — g2 — 0.25, si = S2 = 4, Ai = 
1, A2 = 1, A3 = 6, ki = 2.05, and k2 = 0.05. 


nonlinear waves could be observed in two-mode nonlin¬ 
ear fibers. As an example, we discuss how to observe 
the vector RWs. One can introduce two distinct modes 
to a nonlinear fiber operating in the anomalous GVD 
regime mil, then the ideal condition could be given 
in detail. Firstly, we assume the nonlinear coefficients 
for the two modes are equal, namely, gi = g 2 = 0.25. 
The GVD coefficients are ai = (72 = 1 in the normalized 
units. The spontaneous development of RW seeded from 
some perturbation should be on the continuous waves 
as the ones in [n m m. Experimentally, one can 
use proper picosecond or nanosecond pulses which are 
close to the ideal continuous waves for RW development 
[3^. The amplitudes of the backgrounds should be equal, 
<^1 = <§2 = 1 in the normalized units. One mode laser’s 
frequency k 2 = 0.05, and the other’s is ki = /c2 + 0.5. 
Then the ideal initial optical shape can be given by 
Eq.(19) and (20) with A = Under these 

conditions, the vector RWs could be observed in the non¬ 
linear fiber. The RW distribution is on the retarded time 
dimension, similar to temporal optical solitons. 

On the other hand, it is well known that there are spa¬ 
tial optical solitons in planar waveguide. The similar con¬ 
ditions can be derived directly through coordinates trans¬ 


formation for the coupled NFS in a planar waveguide • 
The RWs in multi-mode planar waveguide could be ob¬ 
served too. Additionally, considering the studies on two- 
components Bose-Einstein condensate(BEC) [H, [sH, [36| , 
we expect that the vector RWs could also be observed 
with the condition. 


IV. DISCUSSION AND CONCLUSION 

In summary, we obtain B-DAB, ABs-ABs, B-DRW, 
AB-ABRW, and vector RWs solutions for the coupled 
model. The corresponding conditions for their emer¬ 
gence are presented explicitly. In particular, when the 
parameter A, which mainly determines the shape of ini¬ 
tial signals, satisfies Eq.(13), there are ideal initial signals 
existing which can evolve to one RW with other types 
nonlinear waves, namely, B-DRW or AB-ABRW, etc. It 
provides possibilities to observe interactions between RW 
and other type nonlinear waves. This is an interesting 
subject for RW studies, since these phenomena are im¬ 
possible to be observed in the single-component nonlinear 
systems. 

Just under the conditions Eq.(16) and (17), which are 
the certain requirements on both backgrounds of the two 
components, the initial condition approaching the ideal 
one Eq.(19) and (20) with A given by Eq.(18), would 
evolve to the vector RWs with no other type nonlinear 
waves. Erom the ideal condition, one can know that the 
two vector RWs can not exist when the two backgrounds 
are relatively static. This is very different from RW in 
the uncoupled system, for which there is no requirement 
on nontrivial background M- Additionally, dark RW, 
predicted in mixed BEC system through numeric stim¬ 
ulation in [25I, is observed analytically here. Moreover, 
one possible way to observe vector RWs in a two-mode 
nonlinear fiber is presented here. 

Recently, the higher-order modulation instability with 
RWs has been observed in nonlinear fiber optics in j3^ . 
There are many possibilities to observe similar phenom¬ 
ena in coupled nonlinear fiber system. We will continue 
to study this subject in the future work. 


Appendix A: The analytic form for B-DRW AND 
AB-ABRW WAVES 

The RW with other type solitons could be derived as 


'tpl 

= V’lO - 

1 

i(A - A)ff>iff>2 

(21) 


i$iP + i$2p + i$3r 

■4)2 

= 4>20 — 

1 

i{\ - A)$i$3 

(22) 


|$i|2 + |$2P + |$3P’ 
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where 


where 


^i[t,z] = ($ 01+^02 + ^ 03 ) 

xExp ^{9i[t,z]+92[t,z]) 

n - 7A/3 - I{2ki - k2)/3 

_ -y^si _ 


^2[t,z] = ( 


^01 [i, 2:] 
^02[t,z] 


^3[t,z] = ( 


T2-/A/3-/(2A:i-A:2)/3 

y</ 9 1^1 

— r2-/A/3-/(2/ci-/c2)/3 

T2-IX/3-I{2ki-k2)/3 


xExp -{92[t, z] - 29i[t, z]) 

_ -V^S2 _ 

n - /A/3 - I{2k2 - ki)/3 

-V^S2 


^03[i,^]) 


^ 0 l[t,z] 

^02[t,z] 


T2- IX/3- I{2k2-ki)/3 

-\/^S2 — T2-I\/3^2k2-k{)/3 

T2- IX/3- I{2k2- ki)/3 


^03[i,x]) 


xExp 


-{Oi[t,z] - 202 [t,z]) 


i{ki - k2) 

"2 2 ’ 

i{k2 - ki) _ 

"2 2 ’ 

+ —(A — ki — k2)‘^Asz‘^ + {A2 + x43)t 

2 

+ -(A -ki- k2){A2 + A^^z + 

2 

+ -(A — ki — k2)Astz Ai A 2 x43, 

= -4i- - 

^1,2 ^1,2 

M{t,z) = Ast'^+ ^{X - ki - k2)'^Asz^ + A2t 
2 

+ 0 (A - ki - k2)A2Z + iA^z 

O 

2 

+ —(A — ki — k2)Astz + Ai. 


Ki 

K 2 

and 


and 

Between the above expressions, the parameter A is given 
^01 = AiEx.p[TitI t^z 2{X — ki — k2)riz/3], by Eq.(18). 

^02 = {A^t + 2A^It2Z + 2/3A^(^X — ki — k2)z + A 2 ) 
xExp[r2t + It^z + 2(A - ki - k2)r2z/3], 

<ho3 = x43Exp[r2t +/r| 2 ) + 2(A -/ci -/c2)t22)/3]. 


Between the above expressions, the parameter A is the 

solution of Eq.(13). Acknowledgments 
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1 

Ki ,2 Kl^ Q{t,z) 
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